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We present a theory of the anisotropy tuned quantum phase transition between spin nematic and 
spin-Peierls phases in S — 1 systems with significant bi-quadratic exchange interactions. Based on 
quantum Monte Carlo studies on finite size systems, it has been proposed that this phase transition 
is second order with new deconfined fractional excitations that are absent in either of the two phases. 
The possibility of a weak first order transition, however, cannot be ruled out. To elucidate the nature 
of the transition, we construct a large- TV SO(3iV) model for this phase transition and find in the 
N — > oo limit that the transition is generically of first order. Furthermore, we find a critical point in 
the one-dimensional (ID) limit, where two transition lines, separating spin-nematic, ferromagnetic, 
and spin-Peierls phases, meet. Our study indicates that the spin-nematic phase is absent in ID, 
while its correlation length diverges at the critical point. Predictions for 23 Na atoms trapped in 
an optical lattice, where the nematic to spin-Peierls quantum phase transition naturally arises, are 
discussed. 



I. INTRODUCTION 

Conventional Landau-Ginzburg theory describes con- 
tinuous phase transitions by fluctuations of an order pa- 
rameter. A nonzero value of the order parameter sig- 
nals spontaneous symmetry breaking and the presence 
of an ordered phase. Moreover, if several distinct or- 
dered ground states are possible, Landau-Ginzburg the- 
ory generically predicts that ordered phases with unre- 
lated broken symmetries are separated by either inter- 
vening phases or a first order phase boundary. In the 
simplest case, this is illustrated by the free energy for 
two independent order parameters <pi and 



t> 2 asi 



r 

f=- 2 ( 



-§(*?- 



Here, a direct second-order transition in the r-g param- 
eter space between phases with either nonzero (p\ or 4>2 
requires fine tuning such that UxU% = u\ 2 is satisfied. 

Recently, a theory of critical phase transitions has 
been proposed that reaches beyond the Landau-Ginzburg 
paradigm and allows direct continuous transitions be- 
tween distinct broken-symmetry phases^ 3 - In this theory, 
the relevant degrees of freedom are described in terms of 
fields that carry fractional quantum numbers and become 
deconfined only at the critical point. In particular, it was 
argued that a direct continuous transition from a valence 
bond solid (VBS) to a Neel phase 2 -! 3 , or a spin-nematic 
phased falls into this class of phase transitions. For this 
reason the search for such deconfined critical phenom- 
ena (DCP) in simple model systems has attracted much 
attention. 

Several systems have been studied up to now. Based 
on quantum Monte Carlo simulations on the two- 
dimensional spin-1/2 Heisenberg model with an addi- 



tional four-spin interaction ( " JQ model" ) the phase tran- 
sition between a Neel ordered state and a VBS appears 
to be consistent with the deconfined critical scenario.^ 
This, however, has been disputed by other numerical 
studies&2Jii that support a weakly first order transition. 
Another candidate system is the iSp(4)-Heisenberg model 
on a square lattice, which describes spin-3/2 cold atom 
systems and has also been conjectured to harbor a di- 
rect second-order transition between a Neel and a VBS 
stated 

The present approach focuses on the anisotropic bilin- 
ear bi-quadratic spin-1 Heisenberg model ("JK model") 
on a square lattice 

H = ^2 [J&i ' Si +X + K(Si ■ S i+X ) ] 

i 

+A^[jS l -S l+y + A-(S l -S l+ ,) 2 ] (2) 

i 

where the exchange integrals for neighboring spins in 
y direction are reduced by an anisotropy parameter A 
(i.e., J — > AJ and A — > AA" on y bonds). For A = 0, 
this Hamiltonian describes decoupled spin chains, while 
for A = 1 the full square lattice symmetry is recov- 
ered. The parameter range of interest is specified by 
A' < J < 0, which is also thought to be the natural range 
for 23 Na atoms in an optical lattice . 12 ' 13 In this range 
the quadratic term favors ferromagnetic order, while the 
quartic term prefers the formation of singlet bonds. This 
competition is solved in two dimensions (A = 1) by 
having a spin-nematic ground state (which breaks spin- 
rotational symmetry but preserves time-reversal symme- 
try (Sj) = 0) and in one dimension (A = 0) by form- 
ing a dimerized ground state (which breaks translational 
invariance) 4 i 12 i 13 i 14 i 15 i 16 When the anisotropy parameter 
is continuously changed between 1 and 0, quantum Monte 
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Carlo studies by Harada et al4 suggest that the system 
undergoes a Landau-forbidden direct second-order tran- 
sition between the nematic and the dimer phases, which 
further motivated a recently developed continuum theory 
for the nematic-dimer phase transition based on DCP. 5 
However, because of significant finite-size effects, they 
were unable to rule out a weak first-order transition or 
the existence of two successive phase transitions^ 

In this paper, we examine the nematic-dimer phase 
transition in anisotropic spin-1 systems using two com- 
plementary approaches. The first approach is based on 
the bond operator formalism introduced by Chubukov 17 
and is particularly suited for studying the dimer phase 
as pairs of neighboring spins are described by common 
bosonic bond operators. Taking the classical limit of 
this approach by neglecting all quantum fluctuations, one 
can already obtain a good overview over the JK model. 
Complementary to the bond operator method, we then 
construct an SO(3./V) model and study its large- N limit 
(N = 1 is the physical limit). This approach explic- 
itly takes the disordering effects of zero-point fluctuations 
into account. For TV — > oo we find that a first-order phase 
boundary separates the spin nematic from gapped spin- 
liquid phases except at special SU(3iV) symmetric points, 
where the transition becomes second order. In the ID 
(A = 0) limit, the existence of the spin-nematic phase has 
been studied extensively but remains elusive . 17 ] 18 i 19 i 20 l n 
this limit, we find a critical point at J = K < where 
spin-nematic, ferromagnetic and spin-Peierls correlation 
functions diverge. The spin nematic phase, however, does 
not exist in a finite parameter range near J = K < at 
A = 0. Finally, predictions for 23 Na atoms in an opti- 
cal lattice, where the nematic to dimer quantum phase 
transition naturally arises, are presented. 



II. SPIN-1 BOND OPERATOR MODEL 

To gain a simple understanding of the anisotropic JK 
model of Eq. © , consider the classical limit of the bond 
operator model following ChubukoAsi 7 -. Grouping the N s 
sites of the square lattice into N s /2 bonds in a columnar 
pattern, we may reformulate the JK model in terms of 
bosonic operators on these bonds that create and annihi- 
late singlets |00) = s^|0), triplets \l,mt) = *m t |0) (where 
m t = 0,±1), and quintuplets \2,m q ) = q^ |0) (where 



m q = 0, ±1,±2). A constraint of one boson per site is 
then necessary to stay within the physical Hilbert space. 
Allowing the bosons to completely condense (as in Bo- 
goliubov theory) then produces a simple phase diagram 
based on the energy of different possible condensates. 

To carry out such program, we must first re-construct 
the spin Hamiltonian of Eq. (|2|) in terms of the above- 
introduced bosons. On a given bond I = (i,i + x), it is 
easier to work with the generators Lg = Si + Sl+x and 
Mi = Si — Si+$ than with the spin operators directly. 
Direct evaluation of the matrix elements of these opera- 



tors in the total spin basis of bond I then tells us how 
to represent them in terms of singlet, triplet and quin- 
tuplet operators. While a little tedious, the net result is 
(dropping the bond labeling £): 



V2(t\t Q + ijt-i) + 2{q\ qi + ql iq - 2 ) 



+ V6(qtqo + qU-i) 

L- = v/2(*S*i + ttifo) + 2(^*22 + ql 2 Q- 
+ V6(qlqi + glx^o) 



(3) 



(4) 



(5) 



\l^(tis + sh ) + y -(<Z^o+4<7o) (6) 
q\h +t\q x +q i _ 1 t- 1 +t t _ 1 <?-i 



M+ = - st< -0 - 2(92*1 - * f _i<7-2) 

- V2( q \t - tU-i) - yf («<S*-i - *t<7o) 



(7) 



M- = --/=(*!« " s f *-i) - 2(4*1 - tl iq - 2 ) (8) 

- V2( q \t - *J*ar-x) - J\ («5*-i - *i<zo) 

To compute the condensate energies for various phases, 
it is then useful to group the bosonic operators into a sin- 
gle vector ip = (s, t\, to, t—i, q 2l q~ 2 ) and to express the 
above operators in the compact form L z = ipa^a g^fft 
etc. Substituting these expressions into the Hamiltonian 
then leaves us with the generic form 

i (W) 

(9) 

where and depend on J/K. That this Hamil- 
tonian consists of only one-body and two-body terms re- 
sults from the one boson per bond constraint. Assuming 
the bosons condense ipg tOL = (tpe,a) + Sipe.on we then ex- 
tract the leading contribution to the condensate energy 
by expanding the ground-state energy to zeroth order in 
powers of 5ij)^ a . 

Four phases are of particular interest: a ferromagnetic 
phase with all spins pointing up, a dimerized spin-Peierls 
phase with a singlet on each bond, a "spin nematic" 
phase in which each site is in the state |1,0) and an an- 
tifcrromagnctic phase with |1,1) (|1,— 1)) on the A (B) 
sublatticc. The bosonic condensates for these idealized 
states are: (q 2 ) = 1 in the ferromagnetic phase; (s) = 1 
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FIG. 1: (Color online) Classical phase diagram. All phases 
are described by a finite condensate of bosonic bond operators 
following Chubukov (Ref. [l7h . Note the existence of a spin- 
nematic phase in the one-dimensional (A = 0) limit. 



in the dimer phase; (go) = \/2/3 and (s) = l/y/3 in the 
spin-nematic phase; and (s) = l/VS, (to) = (—l) v /\/2 
and (go) = 1/V6 in the antiferromagnetic phase (the sign 
here alternates along a column in the columnar dimer 
pattern) . 

Using the above condensates as a guide, we construct 
a simple phase diagram by minimizing the condensate 
energy as a function of (s), (t mt ), and (q mq ) on at most 
two independent neighboring bonds (that is, we explore 
an 18-variational-parameter space). This approach fol- 
lows that of Ref. [l?] in essence and can be thought of 
as a simple two-site clustering method. The resulting 
phase diagram is depicted in Fig. Q] All phases except 
the ferromagnetic phase have a finite dimerization due 
to the explicit translational symmetry breaking of this 
approach. 

It is worth noting the finite existence of a spin-nematic 
phase in the one-dimensional A = limit as shown in Fig. 
[TJ This phase is not expected to survive in the presence 
of quantum fluctuations (finite Stp^a) as pointed out by 
Chubukov-ii. However, the resulting phase diagram after 
these effects are included is not clear from this approach. 
The absence or presence of a (gapped) spin nematic in the 
ID limit has been the focus of extensive numerical stud- 
ies based on density matrix renormalization group, exact 
diagonalization, quantum state transfer, etci 14 ' 19 i 20 ' 21 i 22 
This issue still remains to be settled as diverging corre- 
lation lengths near the ferromagnetic phase limit numer- 
ical approaches. We will show below that in the large-iV 
limit of the JK model the spin nematic phase does not 
occupy a finite parameter region near the critical point 
J/K = 1 for A = 0, while, approaching the critical point, 
spin-nematic correlations certainly diverge. This sug- 
gests that the one-dimensional spin-nematic phase van- 
ishes even for N = 1, although 1/N fluctuations may 
qualitatively change the universality class of the phase 
transition from that of the large- A" limit. 



III. SPIN-1 SCHWINGER BOSON MODEL 

In contrast to the bond operator approach, large-A^ 
methods include the disordering effects of quantum fluc- 
tuations. For the present system, one therefore might ex- 
pect that fluctuations will help restore the translational 
invariance in the dimerized nematic phase and affect the 
nature and the location of the phase boundary between 
the nematic and the dimer phase. Embarking on the 
large-A" route, we introduce bosonic SU(3) spinors ai >a 
(a = 1, 2, 3) to rewrite the spin operators on each site i 



.,23 



Si 



(10) 



In this representation, a\ creates a particle in a state 
whose spin component in the a direction is zero, that 
is, for a given quantization axis, we define aj-jO) = 

(|1, -1)* - |1, l)i)/V2, at 2 |0) = -l)t+\l, l)i)/V2, 

and Q.J 3 1 0) = |l,0)j. Imposing single occupancy on av- 
erage on each site via a Lagrange multiplier fi, one finds 
for the anisotropic JK model (X x = 1, X y = A) 

H= X p[(- J + K ) a i a a i+p, a a iJ1 a i+P,l3 ( n ) 

i 

+Ja\ >a a i+p , a a\ l3 ai + p£ - M^Z ( a \, a a ^ - l)> 



where the summation over the spinor components 
is implied and the exchange integrals are commonly 
parametrized by an angle 6, i.e., J = cos(6*) and K = 
sin(0). The first term in this Hamiltonian possesses uni- 
form SU(3) symmetry, while the second term has stag- 
gered SU(3) symmetry on a bipartite lattice. As a result, 
global SU(3) symmetry is found for special values of J/K, 
namely, for 6 = -3tt/4, -tt/2, tt/4, 7r/2^i 

The order parameter of a spin nematic is a symmetric, 
tracclcss rank- 2 tensor constructed from the spins, i.e., 



(12) 



This tensor can be expressed in terms of the triplons ai >a , 
which when condensed, (ai_ a ) ^ 0, describe long-range 
order. Up to an SU(2) rotation, one then finds that con- 
densation of one of the spinor components (a^ a ) oc x a 
corresponds to nematic order, and condensation of two 
components (a, >a ) oc x a and (a^) oc xp (a 7^ (3), 
where x* a xp is purely complex, indicates ferromagnetic 
order, while additionally staggered expectation values 



with 



i,a{p)) « X* 



C a (/3) 



on one sub-lattice (i e A) and 



on the other (i e B) are associated with 
anti-ferromagnetic order. 

Now, consider generalizing to a large- N SO(3A^) model 
[with SU(3Af) symmetry at = — Stt/A, —tt/2, etc.] and 
expanding in powers of 1/N. To leading order, we obtain 
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a mean-field theory with self-consistent equations 

X P = (jV zZ 4,a a i+P,<x) ( 13 ) 
a=1..3JV 

o*=1..3iV 

where p — x or y. These fields describe short-range cor- 
relations and we choose them to be real to ensure an 
expected time-reversal invariance in the ground state. In 
addition, we restrict the Hilbert space to the SU(3iV) 
representation given by the Young tableau of N columns 
and 1 row by imposing the constraint ^ Q a\ a a iM = N 
(see Ref. [25h . For the following analysis we also allow 
for the uniform condensation of one of the boson species 
dj,! = s/Nx to study the nematic-dimer phase transition. 
The quantity \x\ 2 then describes the nematic "superfluid" 
fraction. 



A. Phase diagram 

Diagonalizing and minimizing the free energy, we ob- 
tain the phase diagram for the JK model in the range 
-0.75tt < 9 < -0.5tt as shown in Fig. EJa). For suffi- 
ciently large A the ground state is always a spin nematic. 
However, as one tunes A towards smaller values the ne- 
matic gives way to disordered spin liquid phases via a 
first order transition except at the tricritical points at 
9 = — 0.5-7T and — 0.757T, where the transition becomes 
a continuous one. This behavior is corroborated by the 
energy gap along the phase boundary [Fig. (2Kb)]. The 
gap vanishes as one approaches the points with enlarged 
SU(3) symmetry at either end of the ncmatic-disordcr 
phase boundary at 9 = — 0.75-7T and — 0.57T. In one di- 
mension [see inset in Fig. [2(b)], the behavior of the gap 
near 9 = — 0.75n can be described by 

A = 6V0/T + O-75exp| m ) (15) 

V y/6/TT + 0.75 J 

where b' and m are positive constants. It is in- 
teresting to note that this form is reminiscent of 
the Berczinski-Kostcrlitz-Thouless-type transition in 
the (l+l)-dimensional sine-Gordon models Equation 
(|15p can be obtained from the constraint equation 
(A -1 a a\ a a,k x . a ) — N by considering the most sin- 
gular contribution in the k x = and 9 = — 0.757T limit 
and by exploiting the fact that the fields approach a com- 
mon value i~j x = Xx = 1/3 near 9 = — 0.757T. Here, A 
denotes the length of the spin chain. 

Our results also show no indication of a gapped ne- 
matic phase near the ferromagnetic phase in the one- 
dimensional regime as was proposed by Chubukov. In 
fact, it is always possible to satisfy the occupancy con- 
straint without a nematic condensate in one dimension. 
However, at the critical point at 9 = — 0.75-7T and A = 0, 



the nematic phase merges from the finite A regime, as 
shown in Fig. [21 so that the nematic correlation function 
diverges at this point. 

Inside the disordered regime we identify three different 
spin-liquid phases [see Fig. [21(a)]: a one-dimensional Z2 
spin-liquid (XxtVx 7^ ®>Xy = Vy = 0)> a one-dimensional 
U(l) spin-liquid (ri x ^ 0,x x = Xy — Vy = 0), and a two- 
dimensional U(l) spin-liquid (rjx,rj y ^ 0,Xx = Xy = 0). 
All three phases are separated from each other by second- 
order phase boundaries that are located at 9 C m — 0.5517T 
and A c « 0.28, respectively. The decoupling of the spin 
chains below A c is likely to be due to the mean-field 
scheme, which overstates the anisotropy in the quasi-one- 
dimensional limit. Moreover, due to the one-dimensional 
nature, we expect only a weak A-dependence of the phase 
boundary between the Z2 and the one-dimensional U(l) 
phases, which could not be resolved. The phase boundary 
between the one-dimensional and the two-dimensional 
U(l) spin-liquids at A c also does not exhibit any depen- 
dence on the JK angle 9. Detailed results along two 
directions in the phase diagram are shown in Fig. [3] 
Last, it is interesting to note that the "unconventional" 
transition point (9 = — 0.557r) examined previously in the 
QMC studies^ lies in very close proximity to the Z2-U(l) 
phase boundary. 

Going beyond mean-field level, saddle-point fluctua- 
tions arc not likely to change the nature of the first- 
order phase boundary. We expect on the other hand that 
Berry phase effects will lead to spontaneous dimerization 
throughout the Xy = Vy = disordered region, which 
lies beyond mean-field level. This situation is similar to 
SU(iV) antiferromagnets, where Berry phase effects arise 
from nontrivial U(l) gauge-field fluctuations and induce 
dimer ordering in ID and 2Di 25 i 27 In particular, following 
such arguments for the SU(3) (J = 0) point in the JK 
model, we expect a columnar dimer ordering for A < Ax 
[see Fig. (2{a)]. We also expect that the one-dimensional 
Z2 spin- liquid phase with Xx,Vx 7^ is unstable towards 
dimer ordering, which can be understood via a mapping 
to the odd Ising gauge theory in ID.— Note that a two- 
dimensional Z2 spin liquid phase is absent on mean-field 
level even at finite A. It is possible that small modifi- 
cations to this model may reveal a two-dimensional Z2 
spin liquid that, depending on the vison fugacity, may 
or may not be stable towards dimer ordering^ This is 
beyond the present study and will be addressed in the 
near future. 



IV. APPLICATION TO OPTICAL LATTICES 

Spin systems with a large bi-quadratic contribution 
are rarely found in solid state systems. However, mod- 
els with higher-order spin interactions can almost per- 
fectly be realized with ultracold spinor atoms in optical 
lattices. Moreover, anisotropy-tuned phase transitions 
between different spin states can easily be induced by 
changing the optical lattice potential in a particular di- 
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FIG. 2: (Color online) (a) Phase diagram of the anisotropic bilinear bi-quadratic spin-1 Heisenberg model in the regime 
— 0.75-7T < 9 < — 0.5tt (J < K < 0) based on a Schwinger boson mean-field approach. The solid line represents a first-order 
phase boundary except at 9 = — 0.5-7T, while the dashed lines represent second-order phase boundaries. The coordinates of the 
tricritical point X separating the nematic, the Neel (9 > — 0.57r), and the disordered phases are #x = — 0.57T and Ax ~ 0.323. (b) 
Energy gap A (in units of f ) along the spin-nematic- disorder phase boundary, which vanishes at 9 = — 0.57T and — 0.757T. The 
inset shows the energy gap for a one-dimensional (A = 0) spin chain of various lengths A near the critical point at 9 = — 0.757T. 
A fit to the largest spin chain (black dashed curve) gives b = ln(6') = —0.945 and m = 0.945 (see main text). 



rcction, hence favoring or disfavoring exchange processes 
between spins along this direction. The spin-nematic- 
dimer transition also conserves magnetization, which is a 
fundamental constraint on optical lattice experiments. 

Considering the lower S = 1 hypcrfine energy mani- 
fold of 23 Na atoms, a large bi-quadratic spin interaction 
naturally arises in an optical square lattice with a single 
average occupancy per sitej 12 i 13 This follows from the 
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FIG. 3: (Color online) Mean field values Tjx,XxjVv>Xv> " su ~ 
perfluid fraction" \x\ 2 , and energy gap A (in units of I) for 
cuts (a) parallel and (b) perpendicular to the 9 axis of the 
phase diagram [Fig|2](a)]. 



two-dimensional spin-1 Bosc-Hubbard model, 

H BH = -* E (U' + J jA') + Y&(' 



(16) 



where b\ a creates a boson on site i with spin a = 0,±1, 

rii = Y^o-H a bi !<7 is the number operator, and Si denotes 
the total spin on site i. The spin-dependent potential U2 
originates from the difference in the scattering lengths for 
the two possible scattering spin channels S = and 2 for 
spin-1 atoms. To leading order in t/Uo, the effective spin 
model then results in Eq. @ (with A = 1), where the ex- 
change integrals are determined by J = —2t 2 /(Uo + U2) 
and K = -2t 2 /(3([/ + U 2 )) - 4t 2 /(3([/ - 2U 2 )). Typi- 
cally, for 23 Na atoms one has U2/U0 « 0.04, which trans- 
lates to J/K « 0.34, or equivalently, 6 « — 0.6047T in the 
effective spin model. In an anisotropic lattice, the hop- 
ping amplitude t becomes directionally dependent, giving 
rise to anisotropic J and K integrals. 

How can one probe the nematic to dimcr phase transi- 
tion? A common way of analyzing optical lattice exper- 
iments is to release all atoms from the trap and to mea- 
sure the column density of the expanding cloud. Within 
this approach, a straightforward way to observe a spin 
nematic to dimer phase transition^ is to apply a weak 
magnetic field (say along the 3-axis) to align the nematic 
hard axis perpendicular to the field prior to the release, 
and to separate different spin components in the expand- 
ing cloud spatially by applying a gradient field after the 
release. A sudden change in the population of the a\ 
state when tuning A then signals the transition to the 
dimerized VBS. 
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Besides this, the interference pattern and the spa- 
tial noise correlations of the expanding cloud also 
provide information about the quantum state in the 
latticej 30 i 31 ' 32 i 33 For instance, in the case of spin-1 atoms 
(with mass m), the equal-time density-density correla- 
tion function for the freely expanding gas after long times 
takes the formal 



G(r,r') = -L £<n a (r)n /3 (r')) 



(17) 



a,/3 



sin 2 (q x L / 2) sin 2 (q y L/ '2) 
L 2 W 4 sm 2 (q x /2)sm 2 (q v /2) 
1 



^e^ R< -^)(S i S J - + (S i S J -)' 



where q = m(r — r')/Ht, L denotes the linear size of 
the square lattice in units of the lattice parameter, and 
W is the width of the expanding Wannier states origi- 
nally centered at lattice sites R^. Here, we have omit- 
ted a delta term from normal ordering and constants of 
order 1/L 2 . The first term stems from the unit occu- 
pancy constraint per site, while the second term contains 
the SU(3) spin structure factor which constitutes the 
JK model exactly at the nematic-ferromagnetic phase 
boundary ( J = K < 9 = — 0.757r). The signature of 
a ferromagnetic condensate therefore is indistinguishable 
from that of a spin-nematic one. Our interest, however, 
lies in the nematic-dimer phase transition. To evaluate 
G(r, r'), we assume complete condensation in the a state 
|0ncm) = Ili a \ cJO) m the s P m nematic phase and obtain 
\' 2 



(s. t s. 



J) 



2 for all i,j. In the dimer phase 

2, 



one gets (SjSj + (SjSj) )dimor = 2 if the spins on sites i 
and j form a singlet and 4/9 otherwise. The correlation 
function then takes the form 



G(r,r') 



+ ■ 



(1 - c)sin 2 (q x L / 2) sin 2 (q y L/ 2) 
L 2 W 4 sm 2 (q x /2)sm 2 (q y /2) 

■(cos(q x ) + COs(q y )) 



W 4 



(18) 



with c = 2 in the nematic and 4/9 in the dimer phase. By 
measuring the density density correlator, both phases can 
be clearly distinguished as shown in Fig. [4] Assuming a 
lattice spacing of 532 nm, a time of flight of 20 ms and 
a detector size of 5 /im, the momentum resolution results 
in Aq/(n/a) ~ 0.015 sufficient to observe the peaks for 
an L w 100 lattice (FWHM/(7r/a) « 0.02, where FWHM 
stands for full width at half maximum). 

Lastly, a recent and interesting proposal for probing 
the spin configuration in an optical lattice is based on 
polarization spectroscopy . 35 i 36 This kind of measurement 
leaves the lattice intact as a signature of the spin state 
is imprinted in the polarization of a probing light beam. 
By analyzing the noise fluctuations in the polarization of 
the outgoing light, the dimerized and the spin-nematic 
phase can be well discriminated. 
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FIG. 4: (Color online) Second-order density correlation func- 
tion for the spin-nematic (blue, dashed curves) and the dimer 
(red, solid curves) phases for different lattice sizes along the 
q y = direction [q = m(r — v')/ht\. The central peak heights 
are roughly given by ~ (1 — c)L 2 with c = 2 in the spin ne- 
matic and 4/9 in the dimer phase. Note the different z-axis 
scalings. 



V. CONCLUSION 

In summary, we have considered the spin nematic to 
dimer phase transition in the anisotropic JK model us- 
ing a bond operator formalism and a large- ./V mean field 
approach. Our analysis generally suggests a first-order 
spin nematic to dimer transition but would not contra- 
dict the possibility of deconfined criticality at J = 
(6 = — 0.57r), where the JK model has an enlarged SU(3) 
symmetry. Our large- TV analysis reveals a critical point 
at J = K < (6 = — 0.757t) where two phase-transition 
lines, separating the spin-liquid, spin-nematic, and fer- 
romagnetic phases, meet. Although the nematic phase 
vanishes in ID, the nematic correlation length diverges 
at the critical point. We have also argued that the ID Z2 
and the ID and 2D U(l) spin-liquids are unstable toward 
dimer ordering while the stability of a 2D Z2 spin liquid 
would depend on the vison fugacity. We do not find a 2D 
Zi spin-liquid at mean-field level but will address a route 
to such a phase in the future. Finally, we have discussed 
various ways of observing a spin-nematic to dimer phase 
transition in optical lattices, where the lattice anisotropy 
can be tuned through laser intensities. Such experiments 
eventually provide concrete ways of studying the phases 
and phase transitions of the JK model. 
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